Structures of Si and Ge nanowires in the sub-nanometer range 
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We report an ab-tnitto investigation of several possible Si and Ge pristine nanowires with diameters 
between 0.5 and 1.2 nm. We considered nanowires based on the diamond structure, high-density 
bulk structures, and fuUerene-like structures. We find that the diamond structure nanowires are 
unstable for diameters smaller than 1 nm, and undergo considerable structural transformations 
towards amorphous-like wires. Such instability is consistent with a continuum model that predicts, 
for both Si and Ge, a stability crossover between diamond and high-density-structure nanowires for 
diameters smaller than 1 nm. For diameters between 0.8 nm and 1 nm, fiUed-fuUerene wires are 
the most stable ones. For even smaller diameters (d ~ 0.5 nm), we find that a simple hexagonal 
structure is particularly stable for both Si and Ge. 

PACS numbers: 61.46.-l-w 
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Semiconductor nanowires with diameters of a few 
nanometers can be grown nowadays by vapor-liquid-solid 
P,|3_, solution-growth 0, or oxide- assisted growth meth- 
ods [J| . These nanowires usually depict a crystalline core 
surrounded by an oxide outer layer. Further removal of 
the oxide layer by acid treatment may lead to hydrogen- 
passivated silicon nanowires as thin as one nanometer Q . 
Pristine (non-passivated) silicon wires with diameters of 
a few nanometers have also been produced from Si vapor 
deposited on graphite . The elongated shape of silicon 
and germanium clusters of up to a few tens of atoms, de- 
termined by mobility measurements 0,0) indicates that 
even thinner pristine structures, with diameters smaller 
than 1 nm, can been produced. 

The growth of such small-diameter structures raises 
the question of the limit of a bulk-like description of 
bonding in these nanowires, since for small enough diam- 
eters the predominance of surface atoms over inner (bulk- 
like) atoms will eventually lead to bonds (and structures) 
distinct from those of the bulk system. In the present 
work, we use first principles calculations to investigate 
several periodic structures of silicon and germanium pris- 
tine nanowires of infinite length, with diameters ranging 
from 0.45 to 1.25 nm. The nanowire structures consid- 
ered are based on the diamond structure, fuUerene-like 
structures, and the high-density bulk structures /?-tin, 
simple cubic (sc), and simple hexagonal (sh). 

Our calculations are performed in the framework of 
the density functional theory Q , within the generalized- 
gradient approximation (GGA) P| for the exchange- 
correlation energy functional, and the soft norm- 
conserving pseudopotentials of TrouUier-Martins 01 in 
the Kleinman-Bylander factorized form We use a 

method 12] in which the one electron wavefunctions are 
expressed as linear combinations of pseudo-atomic nu- 
merical orbitals of finite range. A double-zeta basis set 
is employed, with polarization orbitals included for all 
atoms. For the nanowire calculations, we employ super- 
cells that are periodic along the wire axis, and that are 



wide enough in the perpendicular directions to avoid in- 
teraction between periodic images. All the geometries 
were optimized until residual forces were less than 0.04 
eV/A. Total -energy differences were converged to within 
4 meV/atom with respect to orbital range and k-point 
sampling. 

Most of the nanowire structures we consider in this 
work are derived from crystalline structures. At zero 
pressure and temperature, the structure with the low- 
est energy per atom is the cubic diamond (cd), for both 
Si and Ge, and we use it as the reference structure in 
our calculations. As a first test of the methodology we 
employ, we compute the total energy per atom (Etot), at 
zero pressure, for the following bulk phases: cd, hexag- 
onal diamond (hd), /9-tin, sh, sc, bcc, hep, and fee. In 
table I, we show the total energy per atom of each struc- 
ture, AEtot — Etot — -E't'otj relative to the total energy 
of the cd phase, iS^gj. We observe that AEtot is within 
0.20-0.45 eV/atom for the sc, sh, and /?-tin phases, for 
both Si and Ge. Our AEtot results for the sc, sh, and 
/5-tin phases of Si, and for the /3-tin phase of Ge are in 
good agreement with recent GGA calculations [isl . 
As a further test of the reliability of our calculations, we 
consider the structural transition from diamond to /3-tin, 
that occurs for both Si and Ge under pressure. The ex- 
perimentally measured transition pressure is 117 kBar for 
Si and 106 kBar for Ge J^J. By computing the enthalpy 
H = E +pV from constant pressure calculations, we ob- 
tained theoretical values of 124 kBar and 103 kbar for Si 
and Ge, respectively, both in very good agreement with 
the experimental values. These results also agree with 
previous GGA calculations [3 EE 113 • 

We now address the structure and energetics of several 
stable structures of Si and Ge nanowires with diameters 
between 0.45 and 1.25 nm. Three classes of structures are 
considered. The first one derives from the cubic diamond 
(cd) bulk phase, with the nanowire axis oriented along 
(100) or (110). The latter corresponds to the usual ori- 
entation of observed Si nanowires with diameters between 
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TABLE I; Calculated total energies per atom (AEtot), in 
eV/atom, of selected Si and Ge bulk phases, relative to the 
cubic diamond phase. 





hd 


/3-tin 


sh 


sc 


bcc 


hep 


fee 


Si 


0.096 


0.388 


0.395 


0.442 


0.610 


0.652 


0.673 


Ge 


0.021 


0.240 


0.236 


0.236 


0.296 


0.325 


0.315 



3 and 10 nm 18]. The second class of structures derives 
from fuUerene-like structures 0, . The third class 
of structures derives from the high-density /3-tin, sc, and 
sh bulk phases. Our motivation for the study of the high- 
density structures initially came from the observation, in 
our ab initio calculations, of a structural instability of 
a (6,0) silicon nanotube which spontaneously deformed 
into a simple-cubic structure of considerably low forma- 
tion energy, shown in Fig. ^c). This finding led us to 
consider the /3-tin and the sh high-density phases, which 
also have relatively small formation energies in bulk, as 
shown in table I. In the following, we describe the three 
classes of structures: 

Diamond-structure nanowires - We considered several 
nanowires based on the cd structure with diameters be- 
tween 0.5 nm and 1.2 nm, oriented either along (110) 
or (100). These structures were obtained from the cd 
bulk by defining the wire axis along the indicated crys- 
talline direction, and by including atoms that fall within 
a specified distance from the axis. From this initial geom- 
etry, we removed low-coordinated surface atoms and per- 
formed geometry optimization using the ab initio scheme. 
For both Si and Ge, only the two widest wires, oriented 
along (110), remained diamond-like after geometry opti- 
mization. The geometries of these wires, labeled as cdl 
and cd2, are shown in Figs. (a) and (b) in the case of 
Ge. The corresponding structures for Si are very similar. 
Note that both the cdl and cd2 wires undergo recon- 
struction at the surface but retain a crystalline core at 
the central interstitial channel. The diamond structure 
wires with diameters smaller than 1 nm suffered exten- 
sive reconstructions towards amorphous-like structures. 
Among these, we mention a wire oriented along (100), 
with a diameter of 0.9 nm, which after reconstruction 
became corrugated with pentagonal rings at the surface. 
This wire is shown in Fig. [5] (a). 

Fullerene-like nanowires - We considered two 
fullerene-based geometries proposed in Ref. 0|, namely 
the Si2o cage polymer and the Si24 cage polymer. We la- 
bel those structures as ful3 and ful4, respectively. Based 
on the predicted stability of fiUed-fuUerene-like clusters 
[20|. we also considered variations of ful3 and ful4, la- 
beled as f-ful3 and f-ful4, with the inclusion of two extra 
atoms inside each cage. f-ful4 is shown in Fig. El (b) . Its 
structure is corrugated and presents fivefold rings at the 
surface, resembling the wire shown in Fig. [21 (a). For this 
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(d) /5-tin (e) shl (f) sh4 




FIG. 1: Cross sections of selected Ge nanowire structures 
labeled according to the parent bulk phase. In (a) and (b), 
wires derived from the cubic diamond structure, with axis 
along (110) direction; in (c), a simple cubic wire; in (d), a 
/3 — tin wire with axis along bulk c-direction; in (e) and (f), 
simple hexagonal wires with axis along bulk c-direction. 
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FIG. 2: Side view of corrugated Ge nanowire structures. In 
(a), corrugated wire resulting from a structural instability 
of a (100) cubic diamond nanowire; in (b), a fiUed-fuUerene 
nanowire. Inner atoms are shown as green spheres. 

reason, we classify the latter as fullerene-like, and label 
it as f-ful5. We further investigated filled fullerene-like 
nanowires of smaller diameter, which we label as f-fuU 
and f-ful2. These are based on Sii2 and Siie cages, re- 
spectively, with one additional atom in the center of each 
cage. 

High-density nanowires - The simple cubic (sc) Ge 
nanowire is shown in Fig.^c). It shows very little distor- 
tions relative to the bulk structure. The /3-tin nanowire 
has its axis parallel to the bulk c-axis, passing through 
the center of an interstitial channel. While the initial ge- 
ometry of this wire is somewhat similar to that of the sc 
nanowire, the corresponding relaxed geometries shown in 
Fig. m(c) and (d) are very different, due to the substan- 
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tial relaxation of the former. All four simple hexagonal 
(sh) nanowire structures are oriented along the bulk c 
direction, and retain the crystalline order along the wire 
axis after geometry optimization, regardless of the wire 
radius. Two of the sh structures, shl and sh4, are shown 
in Fig. ^e) and (f), respectively. We also considered an 
empty-hexagon variation of the shl structure in which 
the central atom of the hexagon was removed. We la- 
bel this structure as eh. The relaxed Si nanowires are 
structurally similar to the Ge ones shown in Fig. ^ 
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FIG. 3: Total energies per atom, AEtot (in eV/atom, relative 
to cd bulk energy) of Si and Ge nanowires as a function of 
nanowire diameter. Labeled structures correspond to those 
shown in Figs. 1 and 2. explained in the text. Lines show 
the curves obtained from the continuum model (Eq. 1), with 
dashed (continuum) lines for nanowires based on the cd (sh) 
bulk phase. 



The calculated total energies of the nanowires are 
shown in the Fig. |2| where we plot AEtot as function 
of the nanowire diameter (defined as the diameter of 
the smallest cylinder that contains the nanowire). We 
observe that the formation energies of the high-density 
nanowires (based on sc, sh, and /3-tin phases) and the 
fuUerene-like nanowires, with diameters of the order of 
1 nm, are very close to the energies of the diamond- 
structure nanowires cdl and cd2, of similar diameters. 



with energy differences of ~0.05 eV/atom or less. Note 
that these values are one order of magnitude smaller than 
the energy differences of the corresponding bulk phases 
in Table HJ This suggests that the energetics of wire for- 
mation in this diameter range is strongly affected by the 
surface atoms jU- 

Figure 13 also shows that the amorphous wires derived 
from the instabilities of thin, diamond-structure wires, 
have higher formation energies than those of the high- 
density and the fuUerene-like wires of similar diameter. 
This suggests that amorphous wires could only be pro- 
duced in conditions very far from thermodynamical equi- 
librium. Cage structures like the unfilled fullerene struc- 
tures and the empty-hexagon structures also have very 
high energies as compared to denser structures. The fig- 
ure also shows that, among the nanowires of small diam- 
eters (smaller than 0.7 nm), the shl structure appears 
below and to the left in the energy vs. diameter dia- 
gram, which suggests a high stability for this structure, 
when compared to the other small-diameter geometries. 
In the diameter range between 0.7 and 0.9 nm, fiUed- 
fuUerene-like wires are the most stable ones. 

The above results suggest that the energetics of the 
formation of nanowires are determined by the interplay 
between the energy per atom of a "bulk" part of the wire 
and the corresponding energy of a "surface" part of the 
wire. Motivated by that, we elaborated a simple con- 
tinuum model in terms of bulk and surface energies. Let 
us consider a cylindrical nanowire with radius R, number 
density p (number of atoms per volume) , and total energy 
Enw ■ By decomposing E^^ into a contribution due to the 
bulk-like atoms and a contribution due to the low coordi- 
nated surface atoms, we write Enw = SbpVb+SspVs where 
Eh is the total energy per atom of the bulk-like atoms. Eg 
is the total energy per atom of the surface atoms, and 
Vb {Vs) is the volume occupied by the bulk- like (surface) 
atoms. Vb and Vg are determined by constraining the sur- 
face to a monoatomic layer, such that the surface thick- 
ness becomes The total energy per atom, is 
then written as 



En 



Nat 



+ (Sb - £s) 



[R- 



-1/3)2 



i?2 



(1) 



for R > /O^^/"^, where Nat is the number of atoms of the 
nanowire. For R < p^^l'^ , Snw — £s- 

To test our model, we apply Eq.^to nanowires based 
on the sh and cd bulk phases of Si and Ge. The param- 
eters Eb and p for the bulk phases are obtained from our 
ab-initio results for each bulk phase. The values of Es for 
the sh and cd phases are obtained from best fits to the 
first-principles calculations of the nanowires. In Table ITU 
we show the fitted values of Eg for the cd and sh phases of 
Si and Ge, as well as the difference Ae = Es—Eb for each 
phase and material. The Table shows that Ae is about 
twice as large for the cd phase than for the sh phase 
for both Si and Ge. This means that the energy cost 
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TABLE II: Surface energy relative to the cd bulk energy (sa — 
el'^), in eV/atom, and relative surface energy within a given 
bulk phase (Ae = Ss — et), in eV/atom, for sh and cd-based 
nanowires of Si and Ge. 







(ef - ef) 


(ef - et) 


Si 


0.810 


0.930 


0.415 


Ge 


0.475 


0.530 


0.239 



of a cd surface is much larger than that of a sh surface. 
This probably arises from the fact that surface atoms in 
a cd structure are undercoordinated (coordination three 
or less) , while the surface atoms in a sh structure are still 
highly coordinated, which reduces the energy cost of the 
surface. 

In Fig. 13 we plot Enw given by Eq. 1, relative to the 
total energy per atom of the cd bulk phase, as a func- 
tion of the nanowire diameter for cd-based and sh-based 
nanowires of Si and Ge. The first-principles results for 
the sh phase fall very near Enw, which is a indication 
of the good quality of the model. We also notice that, 
although Es for the cd phase was fitted only to the cdl 
and cd2 wires (the ones that retained the cd structure), 
the Enw curve also passes very near the calculated ener- 
gies of the amorphous wires that originated from the cd 
structure. 

The results of the continuum model seen in Fig. |21 
clearly reproduce the trends of the first principles re- 
sults. In particular, the model correctly describes the 
energy similarity of sh and cd nanowires for diameters 
around 1 nm. The model also predicts an inversion in 
the relative stability of sh-based and cd-based nanowires 
for a diameter Dc around 1.1 nm. This stability inver- 
sion results essentially from the larger for the cd phase 
as compared to the sh phase, for both Si and Ge. It is 
interesting to notice that the structural "transition" of 
the cd nanowires, from crystalline-like (for d > 1 nm) to 
amorphous- like (for d < 1 nm), occurs at diameters very 
near the stability inversion predicted by the model. This 
suggests that both effects are related, as the cd phase 
becomes "metastable" for d < 1 nm. 

Finally, let us comment on the stability of the filled- 
fuUerene-like structures and the shl structure for very 
small diameters. Although Eq. 1 cannot be applied to 
the filled-fullerenc structures (there is not a bulk struc- 
ture associated with them), their calculated energies be- 
have as a decreasing function of the diameter. Such a hy- 
pothetical curve would cross that of Enw for the sh phase 
for diameters near 0.7 nm, as one can see from Fig. O 
This is consistent with the special stability of structure 
shl for very small diameters, and it suggests that the sh 
phase might be the stable one for ultra-thin nanowires. 

In summary, we performed first-principles calculations 
of infinite, periodic wires of Si and Ge with diameters 



between 0.45 and 1.25 nm. We found that the diamond 
structure nanowires are only stable for diameters larger 
than 1 nm; thinner diamond-like wires undergo consider- 
able structural transformations towards amorphous-like 
wires. We propose a continuum model to explain the 
energetics of the nanowires, on the basis of the compe- 
tition between bulk and surface energies. According to 
this model, parametrized from the first-principles calcu- 
lations, high-density nanowires become more stable than 
the diamond structure wires for diameters smaller than 
1 nm. This is consistent with the structural instabilities 
of the diamond-structure nanowires. For diameters be- 
tween 0.8 nm and 1 nm, filled- fuUerene wires are the most 
stable ones. For even smaller diameters (d ~ 0.5 nm), we 
find that the simple hexagonal structure is particularly 
stable for both Si and Ge. 
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